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Abstract 


We have recently presented the complete 0(0^) corrections to the semilep- 
tonic decay width of the b quark at maximal recoil. Here we discuss various 
technical aspects of that calculation and further applications of similar meth¬ 
ods. In particular, we describe an expansion which facilitates the phase space 
integrations and the treatment of the mixed real-virtual corrections, for which 
Taylor expansion does not work and the so-called eikonal expansion must be 
employed. Several terms of the expansion are given for the O(a^) QCD correc¬ 
tions to the differential semileptonic decay width of the 6-quark at maximal 
recoil. We also demonstrate how the light quark loop corrections to the top 
quark decay rate can be obtained using the same methods. We briefly discuss 
the application of these techniques to the calculation of the correction 

to zero recoil sum rules for heavy flavor transitions. 


I. INTRODUCTION 

Precise determination of \Vcb\-i a parameter of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, is an important goal of many experimental studies. The current experimental limit 

i 


lUcfel = 0.036 to 0.046 (90% CL) (1) 

is based on measurements of the beauty hadron decays produced at the T(4S') resonance 
(by ARGUS and by CLEO II) and in Z-boson decays (by the four experiments at LEP). 
In the future large samples of the 6-hadrons collected at R-factories (at SLAG and KEK) 
and at the hadron colliders will increase the statistical accuracy to a few percent level. To 
fully exploit the anticipated experimental improvement, the theoretical description of the 6 
decay must be known with comparable precision. 

There are two methods of extracting the value of |Ucfe|; based on measurements of the 
exclusive decay B D*Ipi and of the inclusive semileptonic decay width of 6-hadrons Pgi. 
These two methods rely on very different theoretical considerations and experimental proce¬ 
dures and complement each other. Their merits and theoretical uncertainties are summarized 
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e.g. in Ref. 0,^. One of the major sources of the theoretical error are the perturbative QCD 
corrections at the two loop level. For the exclusive decays at the zero recoil point these 
corrections have recently been calculated 0,0]. This has significantly improved the accuracy 
of the theoretical prediction for the exclusive method. 

Regarding the inclusive method, recently the O^a^) corrections to differential semilep- 
tonic decay width of the 6-quark at maximal recoil have been calculated |0|. Combined with 
the previously obtained value at zero (minimal) recoil, these results permitted an estimate 
of the complete 0{al) correction to the total semileptonic decay width of the 6-quark. 

In the paper 0] we presented the results of that calculation and discussed its phenomeno¬ 
logical relevance. The purpose of the present paper is a detailed description of the methods 
employed in that calculation. We would like to note that a complete calculation of two-loop 
corrections to a fermion decay width has never been performed before, neither in QCD nor 
in QED (a longstanding example are the two-loop QED corrections to the muon life time). 
Therefore, in the calculation we describe here we had to go beyond the traditional methods 
used in higher-orders calculation. We hope that a description of some technical aspects and 
methods will be of interest for the community. 

Let us hrst mention the difficulties one encounters when trying to compute the fourth 
order corrections to the fermion decays. One of the problems is an appropriate treatment 
of the real radiation of one or two gluons. For the virtual radiative corrections there exists 
a number of algorithms, permitting an efficient, analytical treatment of a large number of 
complicated diagrams, which typically appear in such calculations. On the contrary, no 
similar algorithms were available so far for the treatment of the real radiation. 

The reason why the real radiation at order O^a^) is difficult to evaluate is that the 
particle in the initial state (the decaying 6 quark) carries a color charge and therefore can 
radiate. It is the presence of the massive propagator of the initial quark which makes the 
integrations over the phase space very tough. The kinematical conhguration, where the 
invariant mass of the leptons is equal to zero and the quark in the hnal state is massive is 
the hrst case where a complete analytical evaluation of the real radiation of two gluons in a 
decay of a fermion turns out possible. 

Another potential problem is the treatment of diagrams which represent one-loop virtual 
corrections to a single gluon emission in the 6-quark decay. The virtual corrections in such 
situation are one-loop. Therefore one might naively expect that this case does not require 
any sophisticated investigation. Unfortunately, the integration of the one-loop formulas 
(especially of the boxes) over the three body phase space is difficult. However, it turns out 
possible to express the result for the loop as an expansion which can be easily integrated 
over the phase space. With a systematic algorithm for the expansion, this approach shifts 
the burden of the calculation to the computer. 

The idea which permitted us to calculate the contribution of the real radiation of one 
or two gluons is (qualitatively speaking) the expansion in the velocity of the hnal quark. 
In the limit rric —>■ mi, the charm quark in the hnal state is a slowly moving particle, with 
spatial components of its momentum of the order of — me, much smaller than its mass. 
The momenta of gluons and of leptons (in the case when the invariant mass of the leptons 
is zero) are also of the order of mh — me- It turns out that by a proper choice of the 
phase space variables one can systematically expand the amplitudes and the phase space 
in 6 = (mfe — me)/mb -C 1 (sometimes we also use an equivalent expansion parameter 
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(3 = {ml — ml)/ml = 5(2 — 5); throughout this paper we use as the unit of mass, putting 
ruft = 1). The details of the phase space parameterization will be explained in detail. 

A word of caution is in order here. Our technique proved to be very useful in problems 
where the mass of the quark in the dual state does not differ too much from the mass of 
the quark in the initial state. In this respect an ideal application are semileptonic 5 —> c 
transitions. On the other hand, for such problems as muon or top quark decay, the expansion 
parameter may be close to unity and it is not clear at present if our procedure is of any 
use there. We will, however, show an example where our procedure remains meaningful and 
delivers reliable predictions even in the case when the expansion parameter equals unity. 

The paper is organized as follows. In the next section we discuss how the real radiation 
of two gluons can be computed. Section |TT1| is devoted to the treatment of tensor integrals. 
Then we present a detailed study of one-loop virtual corrections to a single gluon emission 
in the 5-decay. We show that in such diagrams a new type of Feynman integrals appears 
and discuss their evaluation. Section [V] is devoted to the applications of our techniques. 
We present the result for the 0{al) correction to the differential semileptonic decay width 
b ^ c I ui up to 5^^. We also demonstrate how the known light quark corrections to the top 
quark decay width can be obtained with good accuracy if the expansion up to high powers of 
5 is performed, and discuss an application of our technique to corrections to the zero recoil 
sum rules. In the last section we present conclusions. 


II. EMISSION OF TWO GLUONS IN SEMILEPTONIC B DECAYS 

We hrst consider the radiation of two real gluons in a semileptonic 5-decay 

5 —>■ c + e + z/e + (71 + (72 

and discuss how its contribution to the width can be computed. 

Throughout this paper we denote the particles and their four-momenta by the same 
letters, i.e. the momentum of the 5-quark is 5 etc. Moreover, we denote the momentum of 
the lepton pair e -f- z/g as W and often speak about the decay 5 —> c -f W. All formulas in 
this paper apply to the case when the invariant mass of the leptons is zero, i.e. = 0, if 
not stated otherwise. 

Our aim is to construct an expansion in 5 = (m/, — mc)/mb 1. In this case the spatial 
momentum of the c quark is of the order of — me- The momenta of the lepton pair and 
of the gluons are also of the order of mb — me because = gf = gl = 0. The smallness of 
these momenta permits an expansion in terms of 5. 


A. Expansion of the propagators 

To show that such an expansion is indeed possible, we list here all propagators of the 
virtual particles which appear in a calculation of the two gluon emission in the semileptonic 
5-decay: 
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Pi 

P 3 

Pa 

Pe 


{b-giy-ml 

1 


- 2 bgi ’ 


P, = 


{b - g 2 y - ml 


- 2 bg 2 ’ 


{b-gi- g 2 Y - ml - 2 b{gi + 5 ( 2 ) + 2 gig 2 
1 1 „ 1 


ic + gi)‘^-ml 2 cgi 
1 


A = 


1 


{c + gi+g2Y-ml 2c{gi + g2) + 2gig2 


{c + g2Y-ml 2cg2 


{gi + fi' 2 )^ ‘^gigi 


( 2 ) 


All these propagators are shown in the examples in Fig. |l]. If we eliminate the momentnm 
of the c qnark nsing momentnm conservation, c = b — W — gi — g 2 -, 'we can expand the 
propagators P 4 ...PQ in the Taylor series with respect to the “small” momenta gi, g 2 , W. 
Explicitly, we have 


1 ^ ( 2 ^iIT + 2 ^i^ 2 )^- ^ 1 ^ { 2 g 2 W + 2 g,g 2 y 

2 bgi^^ { 2 bgiy ’ ® ‘^bg 2 ( 2 &^ 2 )^' 

1 ^ [ 2(^1 +^ 2 ) 1 ^?' 

2 &( 5 'i + 5 ' 2 ) - ‘2gig2 [ 2 b{gi + 5 ( 2 ) - 2gig2y' 


From these expressions we see that the expansion does not generate any new infrared di¬ 
vergences and, therefore, appears to be permissible. After this expansion only fonr different 
types of denominators remain. We list them for clarity: 

pO _ ^ pO _ ^ 

^~ 2 bgy ^~ 2 bgy 

pO _ __ p _ ^ fp\ 

^ ~ 26 ( 5-1 5-2) - 25-15-2 ’ ^ 2515-2 ■ 

As will become clear from the discnssion below, it is more convenient to use P 3 as one of 
the basic propagators than to expand it any further. 

Let us emphasize at this point the advantages of the above expansion. As can be seen 
from the above expressions, the four-momentum of lepton pair W and the momentum of 
the c-quark do not appear in the denominator. This immediately implies that the integral 
over the phase space of W and c factorizes. Therefore, the expansion suggests a simple way 
how the rather nontrivial integration over the phase space of four particles in the hnal state 
can be reduced to a more familiar case of the three body phase space. This reduction is 
described in the next section. 


B. Phase space integration for the emission of two gluons 

After we have checked that the propagators of all virtual particles can be expanded, the 
same must be demonstrated for the phase space element. This is done in the present section. 

Considering the semileptonic decay of the 6-quark and integrating hrst over the lepton 
phase space, one obtains the following phase space integration element (we consider only 
the point where the invariant mass of the lepton pair IF^ is zero): 

di ?4 = [dc][dW][dgi][dg 2 ]b^{b - c - 51 - 52 - W), 
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where a shorthand notation is 

= D = ^- 2 e. 

It is convenient to introduce an auxiliary vector H which equals to the sum of the momenta 
of the c-quark and W boson, H = c + W. With this notation the phase space integral 
becomes 

= J dH^ J[dH][dgi][dg 2 ] 6 ^{b - H - gi - g 2 ) j [dW][dc]5^{H -W-c). 


The integration over the (IT, c) phase space gives (for now we assume the integrand contains 
no IT in the numerator; the tensor integrals with IT dependence will be analyzed in section 


m 


jdWc 

1 X 2 


j[dW][dc]5^{H -W -c) 


VLd-i {H^ - 

2D-1 ('^2)D/2-l 


where 




2 7r3/2-^ 

r(3/2 - 6) 


is the volume of a D — 1 dimensional sphere of a unit radius. 

Having performed the integration over the (IT, c) phase space we are left with the three 
particle phase space of one massive (if) and two massless {gi, ( 72 ) particles and the integration 
over the square of the momentum if. 

We use the variables (with ml = 1) 


= y, {91 + 92 ? = U, (if + g 2 ? = z, 


(5) 


and get the following expression for the four particle phase space []1 
2^ j B{l-e 1 - e) / 

with Um = {I — z){z — y)/z and the following limits of integrations 

ml < y < 1, y < z < 1, 0 < u < Um- 

In order to expand the phase space element in powers of /? = {ml ■ 
change the variables in the above expression 


{y - ml) 


2\l-2e 


y 


1-e ’ 


( 6 ) 


(7) 


ml)/ml it is useful to 


y=l-(3xi, Z=l-(3xiX2, U = UmX3, 

_ f3‘^x\x2{l - X 2 ) 

1 - (3xiX2 

The integration limits for all Xi are the same: 

0 < Xj < 1 


(8) 
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Expressing the phase-space element through these variables and restoring proper dimen¬ 
sionality we End: 


di?4 = Rl 


/■l 

1 1 

1 U.XU.X2hJ.X3 U *^2? *^3/ 

(1 - (3xiX2)^~^ (1 - (3xi)^-\ 


where 


J{xi,X 2 ,X‘i) = x\ 


\l-2e l-2e 

Jbr) 


(1 -X 2 ) 


l- 2 e -e 

Jjn 


(1 


Rl = 


[m. 


2\2-3e 




5(1 - e, 1-e) V 2^-1 


(9) 


Our hnal aim is to expand di ?4 in the limit of The above expression is well 

suited for this purpose. Expanding it in powers of /3 = {ml — ml)/ml gives rise to simple 
integrals which can be expressed by Beta functions: 


dxi da;2 <lxzJ{xi,X2.,xz)xix'{^x^ 

= 5(n + 4 - 4e, 2 - 2e)5(m + 2 - 2e, 2 - 2e)B{k + 1 - e, 1 - e). 


( 10 ) 


C. Basic propagators 


The basic propagators P°, P 2 > R^) P 7 , expressed in terms of the variables yield 

simple expressions which can be integrated over the phase space. We hnd 


P/ 

P 3 


1 _ 1 
2bgi (3 xiX2 ’ 
1 


2bg2 

1 


2 &(fi'i + 5'2) - 25-15-2 I3xi 


1 - I3XiX2 

/3a;i(l - X2)[l - f3xiX2{l - xa)]’ 

3 _ 1 _ 1 - (3X1X2 

^ 25 ( 15(2 ( 3 ‘^XIX 2{1 - X 2 )X:i' 


( 11 ) 


Clearly, the expansion of these quantities in terms of small (3 is readily performed. The 
resulting integrals are very similar to the integrals in eq. (|^) and can be expressed in terms 
of the Beta functions. 

The same is true for all scalar products of the four momenta which enter the calculation. 
Therefore, the above discussion demonstrates that the part of the semileptonic decay width 
of the 6 -quark containing radiation of two gluons can be calculated by expanding the matrix 
element and the four particle phase space in powers of (3. What is still missing in our 
discussion is the study of the integrals over the (IT, c) phase space when the numerator 
depends on W. We discuss this issue in the next section. 


III. TENSOR INTEGRALS 

The expansion of the propagators in terms of the momenta IT, gi, g 2 , may result in high 
powers of the scalar products (Wgi)°'(Wg 2 )^, which have to be integrated over the (IT, c) 
phase space. In this section we present efficient methods for such integrations. The results 
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of this section are not restricted to the case = 0 and are likely to be useful for solving 
other problems. 

The main object we are going to discuss is the following integral over the (IT, c) phase 
space: 


I(a,b) = I [dW][dc]6^{H -W- c)(Wgi)yWg 2 f 

= ((WginWg2)^}H=w+c ( 12 ) 

for arbitrary a, b > 0. This integral can be rewritten as 

I(a,b) = T^^a+b 

= {W^^....W^^)h=w+c- (13) 

The important property of the tensor T is that it depends on a single external vector only 
{H). This property simplihes the calculation of the necessary integrals. We discuss below 
two approaches to such calculation; their relative merits depend on the capabilities of the 
symbolic manipulation languages, if the algorithm is to be implemented using a computer. 


A. Method 1 

The hrst approach is more general and potentially more efficient. To fully exploit the 
dependence on the single external momentum we can extract the component of the vector 
W which is parallel to it. For this purpose we write: 

WH 

W = W^ + c^H, W^H = 0, c, = -— (14) 

The value of Ci is then hxed, because in the general case IT^ 7 ^ 0 we have 

WH =^{H^ + W^ -ml). (15) 

Therefore, if we use the above substitution for W in the tensor integral we obtain a 
sum of tensor integrals over transverse components of the IT (with respect to H): 

Tr = (fTf Wt)H=w+c. (16) 

Evidently, the tensor structure of this integral is trivial: since Wj^H = 0, Tj)’"' must be 
proportional to the absolutely symmetric tensor of the {D — 1) dimensional space. The 
“building block” is the metric tensor of {D — 1) dimensional space: g^^^, — H^Hy/H'^. We 
note that only tensors of even rank contribute. The generation of the absolutely symmetric 
tensor of H —1 dimensional space can be easily encoded in a symbolic manipulation program. 
Hence, the algorithm described above permits an efficient treatment of the tensor integrals 
which appear in the problem at hand. 
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B. Method 2 


Using the same idea as in the previous section we can advance slightly further with the 
analytical calculation of the integral I{a,b). 

After writing W = W± + ciH and noticing that giH and g 2 H are constants in the (W, c) 
phase space, we conclude that the non-trivial integrals to be computed are 

I±{a,b) = {{W±giy{W±g 2 )^)H=w+c = {{W±gi±)°‘{W±g 2 ±)'^) h=w+c, (17) 

with gi± (i=l,2) dehned in analogy to W±. 

The integral I± is a Lorentz-scalar and we can choose an arbitrary frame for its calcu¬ 
lation. For instance, we consider H to have only the zeroth component and put the ; 2 -axis 
of the D — 1 dimensional space along the vector gi^. Then, rewriting g 2 ± as 

g2± = g2a + C2gi±, = 0 , C 2 = (18) 

we hnd that I± is expressed as a sum of integrals of the from: 

h{a, b) = {{W^gi 2 .T{W^g 2 aY)H=w+c. (19) 


This integral is non-zero only if both a and b are even integers. 

Since g 2 agi± = 0, we can put the x-axis of the D — 1 dimensional space along the vector 
g 2 a- Then we have: 


Ii{a,b) 

h{a,b) 


^ W^g,^ y/ W^g2a V 

J ^D-i [\w^\\gi±\) [\w^\\g 2 a\) ■ 


( 20 ) 


The integration element is 


D-2 


dhl/)-! = K[ sin^ ^ ^ 6i d6i, 0 < < tt {i ^ D — 2), 0 < 9d_2 < 27r 


i=l 


If we choose 2 : and x axes as described above, we get 

W±gi_L 


\w^\\gi±\ 

W±g2a 

\w^\\g: 


= — cos 6 ^ 1 , 

= — sin 9i cos 6 * 2 . 


2a I 


( 21 ) 


( 22 ) 


The minus sign in the above equations arises because the transverse vectors here are space¬ 
like. In calculating 12 ( 0 , 6 ) the following formula is useful: 



1 -\- cn 1 -|- / 3 \ 

2 ’ 2 ) ■ 


(23) 


This formula is valid for even and positive (3. Using it we get the result for the integral 
h{a,b): 



(24) 


h{a,b) 


(_l)a+bF 1 

(i- 


|F| 


|F| 


ttF I 

f CL-\-b-\-3- 
{ 2 

■“J 

1 


In order to use these results for the further integration over the phase space {H,gi,g 2 ) it is 
necessary to express the modulus of the vectors W±, gi±, g 2 a in terms of the variables {xj} 
introduced in the previous section. This is a cumbersome but straightforward task and we 
do not discuss it here in any detail. 

The steps described above allow an easy implementation in any symbolic manipulation 
program. They provide an efficient and uniform treatment of the tensor integrals which 
appear in this calculation. 


IV. ONE LOOP VIRTUAL CORRECTION AND ONE REAL GLUON 

Another source of the 0{a‘^) corrections are the one-loop virtual corrections accompanied 
by a radiation of an additional real gluon. Again, we would like to expand those diagrams in 
powers of 5 = {mb — me)/mb. The purpose of this section is to provide a suitable algorithm. 

There is a principal difference from the radiation of two real gluons in this part of the 
corrections. Namely, it is not sufficient here to perform a Taylor expansion of the relevant 
diagrams. 

To see this, let us analyze a simple example. We consider a diagram where the &-quark 
radiates a gluon and then a one-loop self-energy correction is inserted on the b quark line 
(see Fig. Taking the limit mb —>■ me in such situation implies that the momentum of the 
real gluon goes to zero. Therefore the self-energy diagram should be calculated close to the 
mass shell of the 6-quark. It is well-known that in this limit the self-energy diagram has 
an on-shell logarithmic singularity of the form 

log(M^ — p^) ~ log(M^ — m^) ~ log((5). 

Since this diagram contains log(5) it is clear that it cannot be calculated as a Taylor 
expansion with respect to 6 . Therefore, a more sophisticated procedure is needed and the 
method of eikonal expansions [|],^ is used in this case. We illustrate this approach in the 
following example. 


A. An example of the eikonal expansion 

In this subsection we consider a simple example of a one-loop integral, where a result 
obtained from the eikonal expansion can be compared with an exact formula. We calculate 
the scalar self-energy diagram: 

B(P^ = J ^2Tr)D k^[{k + p)^ - m^y 

Performing Feynman parameterization, integrating over the loop momenta, and expanding 
in p = (m^ — p^)/p^ we get (using p^ = 1 here): 
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B^p^, m^) 


(26) 


^r(l + e) 

(47r)-f^/2 


'1 , p2 3 4 

- + 2-p + p\np- —+ -- — + ... 
e 2 6 12 


The appearance of the log(p) term in this resnlt signals that the Taylor expansion of the 
integrand will be insnfficient. In order to get the correct resnlt we add the eikonal expansion. 
For this pnrpose we write 


B{p^, irf) = B, + B,ik 


(27) 


The Taylor expansion term Bt is obtained by expanding the integrand in (|23| ) in a Taylor 
series in p. This yields a set of one-loop on-shell integrals which are well known. The validity 
of snch an expansion is determined by the condition k‘^ + 2 pk 3> pp^. If this is not the case, 
the Taylor expansion breaks down. This breakdown gives rise to an infrared singnlarity 
which occnrs in the region where pp^ ~ 2pk and -C pp^. This spnrions divergence can be 
cancelled by adding an expansion of the integrand in k"^ (for a more detailed discnssion see 



One shonld note in addition that if any power of appears in the nnmerator we get an 
integral of the form: 


lin) = 

J {2 tt)^ 2 pk — p^p^ 


n > 0. 


An important point is that snch integrals are zero within the dimensional regnlarization 
framework. One simple reason for this is that the integral over transverse components of k 
is scaleless. 

Hence, the only term in the eikonal expansion we shonld consider is the one with /c^ 
neglected in the second propagator. The eikonal integral therefore reads: 


Heik — 


d^k 1 

(27r)^ k‘^{2pk — p‘^p) 


(28) 


It is evident that this integral scales with p as Heik cc It is this dependence on e which 

gives log(p) in the hnal resnlt. 

Consider i?eik in the rest frame of p. Performing the integration over the time component 
of the loop momenta first, one gets a simple integral representation over the transverse 
components of the loop momenta which can be easily calcnlated. More sophisticated eikonal 
integrals are considered in the next section, where all necessary formnlas can be fonnd. 


B. Eikonal integrals: preliminaries 

Here we provide a set of formnlas for dealing with integrals appearing in the eikonal 
expansions. As is clear from the example presented above, the region of integration in which 
we are interested here is characterized by the valnes of the virtnal momentnm of the order 
of the momentnm of the real glnon emitted in the decay. 

We denote the loop momentnm by k and the momentnm of the real glnon by p; explicitly, 
the eikonal expansion of the propagator 

{k + h — gY — ml P + 2 k{b — g) — 2bg 
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is given by 


P{k)eik = 


1 

2kb — 2bg 


OO 


E 


/ 2kg - e 
\2kb — 2bg 


j 


The basic eikonal integral in the one-loop correction then readsj^ 


I^{n, m, a, b) 


d^k {kWY{kgf 

{2ti)^ (}P + i6){2kQ ± 2Qg + i6Y{2kQ + 


(29) 


In the actual calculation Q is always a time-like vector which for different diagrams can 
be 6, b- g, c + g. 

As we pointed out above, the integrals with any power of in the numerator vanish. 
There are several ways to see this. First, the integrals over transverse (with respect to Q) 
components of the loop momenta are scaleless; also if there is no in the denominator, the 
poles of the integrand are located only on one side of the k^ integration axis. Also, if n < 0, 
the integral in eq. (|^) is zero because it is scaleless. 

Let us describe the most economic way to calculate /^. We choose the Lorentz frame in 
which Q = (Q^, 0,..., 0). The integration over the time-like component of the vector k can 
be performed using the residues. If the integration contour is closed in the upper half-plane, 
only that pole in contributes for which fco = —k±. 

After fco) fhs integration over the D — 1 dimensional space should be performed. To 
discuss the integration over fc_i_ we first note that after integrating over k^ the denominator 
of the integral /^(n, m, a, b) depends on k\ only. Therefore, all tensor integrals which appear 
due to the numerator structure are simply related to the D — 1 absolutely symmetric tensor 
of the corresponding rank. Hence, the tensor integrals here are similar to the ones discussed 
in Section B 

Finally, we are left with the radial integration over the modulus of k^. The following 
formulas are useful: 


I (n, m, mi) = J 


d^k 


(kl) 


2 \mi /2 


/■^(n, m, mi) = J 


{27i)^ -|- i6){2kQ — 2Qg + i6)^{2kQ -|- iS)'^ 

-— -j^B{2 — m + mi — 2e,n + m — mi — 2 + 2e) 

(27r)'^ ^ 

d^k 

(27r)'^ -|- i6){2kQ + 2Qg + i6)'^{2kQ + 


iVt D— 


h(Q9) 


2—m+mi —n—2e ^^2^ — 1—mi/2+e 


(27r)^- 

('_^)n+m+mi+l2-m-n-l ^(2 _ ^ _ 2e) 


X 


X 


(n — 1)! F(3 — n — m + mi — 2e) 

^{2 — m + mi — 2e) — V’(—1 + m — mi + 2e) , 


(30) 


^In the Feynman diagrams with self-gluon couplings the situation is more complicated. Their 
treatment is described in subsection IV Q. 
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with 


i){z) = 


1 dT{z) 
r(;2) dz 


In the above formulas IId-i is the volume of the D — 1 dimensional sphere of unit radius. 

There is one subtle point concerning the expressions presented above. Namely, it is 
easy to see that the integrals have an additional singularity after the integration over ko 
has been performed. This singularity corresponds to the appearance of the imaginary part 
in some of Feynman diagrams which contribute to the result. Below we explain how this 
singularity was treated. 

In the singular integrals the integration over the modulus of the transverse momenta is 


Jo {k±-l-i5)^' 

The singularity is located at the point k± = 1. To make this integral meaningful we rewrite 
it as 

1 / d k^ I 


(n — 1)! yda; 


0 


Rescaling k± ^ x k± we get: 


(n — 1)! \ dx 




For the derivative we use 


x^ = 


F(1 + a) 
T{2 + a — n) 


Finally, after performing the remaining integration over kx_ using 


Re / dx-- = ihil + a) — ihi—a), 

Jo x-l + iS ^ ^ 

we arrive at the formula for quoted above. 


C. Eikonal integrals for the graphs with the gluon self-coupling 

The graphs with a triple gluon coupling lead to the most difficult eikonal integrals. 
The difficulty originates from the fact that in such graphs one has two massless (gluon) 
propagators. A typical integral reads 


(27r)'^ {k — gY {{k + b — gY — + kY — rnl) 
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The Taylor expansion of such integral is easy. However, in the eikonal expansion one 
cannot expand massless propagators. Therefore, the master integral in this case is 

_ 1 _ 

J ( 27 r)^ F {k - g)\2kQ - 2QgY{2kQ)-- ’ 

where Q is again a “large” time-like vector, as described in the previous section. 

Note also, that in this case the terms with or {k — gY in the numerator do contribute 
to the result. 

Let us describe in detail how such integrals can be calculated. For this purpose we first 
combine the massless denominators: 

1 /■! dx /■! dx 

{k — gY Jo (/c^ — 2kgxY Jo {{k — gxYY 

Next, we decompose the product of two eikonal propagators 

1 

{2kQ - 2QgY{2Qk)^ 


using partial fractions. The integral I (m, n) becomes now a sum of the integrals of the 
following type: 



d^k 1 

(27r)'^ [{k — gxYY{2kQ — 2gQY 

d^k 1 

(27r)'^ [{k — gxYY{2kQY 


In both integrals we shift the integration momentum k ^ k + xg and get 


(32) 


/i(n) = f 

Jo 

h{n) = f 
Jo 


dx 


dx 


d^k 1 

(27r)^ (F)2[2fcg - 2gQ(l - x)Y' 
d^k 1 

(27r)^ {k‘^Y{2kQ + 2gQxY 


(33) 


Now, in the hrst integral one makes the change of variables a; —*• 1 — x and then rescales the 
integration momenta k ^ xk m both Ji and l2- The integration over x factorizes and can 
be done using 


dx x'^ 


1 


JO a -f 1 

The integrals over k differ from the integral in the previous section by the presence of 
two powers of in the denominator. In this case it is useful to perform a transformation 
(similar to the technique of integration by parts f^), which reduces the integrals with the 
second power of in the denominator to the integrals with the hrst power only. 

To write down this transformation we introduce the notation 



d^k 1 

(27r)^ {k‘^Y{2kQ ± 2QgY{2kQ)^' 
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Then 


/^(2, n, m) = — 2{m + n,m + 2) — 2 n + l,m + 1) 


Using this formula we end up with integrals with the hrst power of k‘^ in the denominator, 
for which the formulas of the previous section are applicable. 


D. Integration over the single gluon phase space 

In the treatment of the radiation of one gluon with 0 (a‘^) accuracy, the last step one 
has to do is to perform the integration over the phase space of the real gluon. The eikonal 
expansion, described above in detail, gives the virtual corrections in terms of powers and 
logarithms of the small parameters. This form of the intermediate result simplihes the hnal 
integration over the phase space of the decay products. This would not be the case if we 
had an exact result for the loop: that would certainly contain dilogarithms of complicated 
arguments. 

The phase space integration with a single gluon is in its general structure similar to the 
case of two real gluons, discussed in section m For completeness, we give a short account 
of necessary formulas. 

In the present case the phase space element is: 


with 



(b-c-g-W), 


(1 - px) 


l-e 


(34) 


H = W + c, = 1-Px, 

= 1 - bg = Hg = —. 

The master integral in this case reads: 

[ dxJ{x)x^ = B{n + 2 — 2e, 2 — 2e). 
Jo 


(36) 


(36) 


The results of the eikonal integrals require a slight modihcation of the master integral. As 
we saw in section [IVB| , those results contain {hg)~‘^'^ which results in an extra factor 
In this case the master integral is 


dxJ{x)x^-^^ = 5(n + 2 - 4e, 2 - 2e). 
Jo 


(37) 
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V. APPLICATIONS 


In this section we discuss some applications of the techniques described above. 

First, we give a complete result for the correction to the differential semileptonic 

decay width of the 6-quark at the maximal recoil point. We next re-derive the BLM 0(ag) 
corrections to the quark top decay width into massless hF-boson and b quark with the help 
of the expansions presented above. As will be explained below, the expansion parameter 
in this case is unity. It not obvious that the techniques described above can be of any use 
there. Therefore it is interesting to check how the procedure works in such extreme limit. 
Finally, we comment on the applications to zero recoil sum rules. 


A. Two loop QCD correction to semileptonic ft^decay at maximal recoil 


The measurement of the inclusive semileptonic decay width permits a determination of 
the CKM matrix parameter \Vch\ with a small theoretical uncertainty [Q,^. The magnitude 
of the perturbative corrections to this quantity has been subject of discussions in the 

recent literature. We addressed this problem in ref. where the exact calculation of the 
0(ag) corrections at maximal recoil was used to estimate the total 4th order correction to 
the b quark semileptonic width. Here we present more complete results of that calculation. 

We consider semileptonic decay of the 6-quark 6 — dui. The momentum carried away 
by leptons is denoted by q. We write differential semileptonic decay width of the decay 
6 —>• dui at = 0 as 


dF,i 


dg^ 


q2=0 


ro 


Aeorn H- Cp^i + 

TT 



Cp^2 


(38) 


where Fq = ^3 |WbP- ABom = (1 — and Ai is also known in a closed ana¬ 
lytical form |IT| , p!^ . A 2 is the new result which gives correction to the differential 

semileptonic decay width of the 6-quark at the point of zero invariant mass of leptons. 

For the purpose of the presentation we write the 4th order correction as 


A 2 — 6^ [{Cp — Ca/2)Af -|- CaAa + TpNpAp -|- T^Ah] ■ (39) 


In the above equation 6 = (1 — rric/mb) is the expansion parameter. Ap describes the 
contribution of the massive 6 and c quark loops. 

For the SU{3) group the color factors are Ca = 3, Cp = 4/3, Tp = 1/2. Np is the 
number of light quark flavors whose masses were neglected. 

In ref. [|] we have presented analytical results for the functions Ap^a,l,h up to 6^, though 
actual calculations have been done up to 6®. Our complete formulas for these functions are 
presented in Appendix ^ 


B. BLM O(Poag) corrections to the decay width of the top quark 

As another application of the above techniques we consider the decay of the top quark 
t — Wb. It is well known that this decay width (at least at the Born and the one-loop 
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level) can be well approximated by neglecting the masses of the hh-boson and the 6-quark. 
In such case, since = 0, the techniques presented in the main part of this paper can be 
readily applied. In particular, the formula ( p8D can be rewritten in such form that it gives 
the correction to the two-body decay width t —>• QW (where Q is a heavy quark): 


Tit ^ QW) = To 


ABorn H- -Cp^l + ( — ) Cf^2 

TT V VT 


(40) 


The Aj functions here are the same as in eq. (pSf) , with mi, and replaced by mt and mg. 
The problem, however, is that the procedures described above were based on the expansion 
of the rate in the mass difference of the hnal and initial quarks. In the case of the top 
quark decay with a b quark in the hnal state this means that the expansion parameter 
S = (mt — mi,)/mt is close to 1 for realistic values of and mi,. 

Little is known about convergence properties of the series described in this paper. For¬ 
tunately, part of the Oia"^) correction to the top quark width is known exactly; it is the 


contribution of the massless quarks, calculated for mw = mb = 0 in 00 (the relevant 
diagrams are shown if Fig. |^). We can use that limiting case to check if we can reproduce it 
with our techniques. 

The exact formula for the massless quark correction reads (in the MS-scheme and for 
one generation of massless fermions) 


F^"™(f ^ Wb) 


CfTr 
TT ) 16 


Gf 

~ 8^271 


28 

2C(3) + y C(2) 



(41) 


It should be mentioned that the diagrams with real or virtual massless fermions represent 
the simplest case for the described algorithms. The reason is their simple planar topology, 
which allows the computer programs to work very fast. For this particular type of diagrams 
we can expand the width up to a very high power in 6. 

We write the result of this expansion as 


Ff"™(t ^ bW) 


(^)'c^T^Fo F(<5), 


where the function FiS) is given as a series of powers and logs of 5; we have computed this 
expansion up to terms of the order 5^^. The hrst few terms of that expansion are given by 
the function in the appendix. 

The numerical value of the exact result is 

F(l) = i(2C(3) + fc(2)-5) = 0.3574.... 

The values of the approximate result Fib) for <5 = 1 are, for three numbers N of the summed 
terms. 


N = 21, 

F(l) = 0.3266, 


N = 15, 

F(l) = 0.3210, 


N = 11, 

F(l) = 0.3176. 

(42) 
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Comparing these numbers with the exact result we see that already the hrst 11 terms of 
the expansion give the accuracy of about 10 %. Unfortunately, the accuracy does not improve 
signihcantly with the growing number of terms in the expansion. This slow convergence of 
the series is caused by the logarithms of the small {b quark) mass. Although these terms 
are suppressed by powers of that mass (and should give no contribution to the hnal result if 
mb = 0), they spoil the convergence of the expansion. It would be very useful (and we think 
it is possible) to hnd a systematic way of eliminating those parts of the integrands which 
give rise to these logs. 

On the other hand, from the perspective of practical applications (like top physics at the 
Next Linear Collider) it would be sufficient to know the two-loop QCD correction to the top 
quark decay width even with 10 percent accuracy. It is tempting to apply our algorithms 
to this problem and extrapolate the result to the point h = 1; on the other hand, at present 
we do not know any reliable method for estimating the uncertainty of such result. 


C. Zero recoil sum rules for 6 — > c transition with 0{ag) accuracy 


Another useful application of the techniques described above is connected with the so- 
called zero recoil (ZR) sum rules |^. In this subsection we would like to discuss this point. 
A more detailed discussion of ZR sum rules with the 0(a‘^) corrections will be given in a 
future publication 
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The ZR sum rules are important for the estimates of the zero recoil transition form factors 
such as Fb^d* for B —> D* transition or Fb^d for B ^ D transition. In turn, the formfactor 
Fb is a crucial ingredient for \Vcb\ determination from the exclusive semileptonic B 
D*li'i decays. 

The sum rules for exclusive heavy-to-heavy flavor transitions are based on the operator 
product expansion (OPE) of the hadronic amplitudes in terms of the inverse quark masses. 
The zeroth order term in this expansion is the parton model where free quarks are substituted 
for real hadrons both in initial and final state. Here we disregard the non-perturbative 
corrections and discuss how the perturbative corrections (P(a^) to the sum rules can be 
evaluated. 

We consider a transition of a 6 quark at rest to a c quark and massless partons which 
occurs under the influence of the external current J^. The momentum carried away by the 
external current is g = (go, 0), where 


% = e < p <C mb, me- 

The quantity which is of primary importance for the perturbative corrections to ZR sum 
rules can be schematically written as: 

Jdewj{e), wj{e) = rj{b ^ c +X \ qo). 

0 

In this equation we have shown explicitly the dependence of the transition rate rj(6 —>■ c | go) 
both on the external momenta transfer go and on the Lorentz structure of the current. 
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If e = 0, the transition is elastic; the hnal state is a single c-quark at rest. The 0(a 


corrections in this case reduce to the renormalization of the external current For vector 


and axial currents they were calculated in 

On the other hand, if e 7 ^ 0 the c quark starts moving and can radiate. The second order 
corrections 0(as) were calculated in 0 . 00 - Aiming at the fourth order, i.e. at 0(ag) 
accuracy, one has to consider the hnal state with the c-quark and two real gluons or light 
quarks and also the 0(as) correction to the single gluon emission in b —>■ c transition. 

Due to the hierarchy of scales, e < /i <C mi,, me, the hnal c-quark is moving slowly by 
dehnition; this is therefore a nice place where the techniques described in this paper can be 
applied. In particular, an algorithm for performing eikonal expansions with the subsequent 
integration over the phase space appears to be very useful here. This topic will be discussed 

ra¬ 


in detail in ref. 


VI. CONCLUSION 

In this paper we have described techniques used in the calculation of the 0{al) correc¬ 
tions to the semileptonic decay width h ^ clui at maximal recoil [^]. 

The technical tool we used for that calculation is an expansion of the decay rate in powers 
and logarithms of the mass difference between the initial and hnal quarks. We presented 
a detailed discussion of the algorithm, which enables one to construct such expansion. We 
treated virtual corrections, emission of one gluon, and emission of two gluons separately. 
Therefore, these algorithms can be used for the analyses of less inclusive quantities than the 
total decay rate, at least in principle. 

In case of two-loop virtual corrections and emission of two real gluons, the expansion 
in 5 = (mb — me)/mb is a Taylor expansion. In the case of one-loop corrections to the 
amplitude of single gluon emission in &-decay, the Taylor expansion is insufficient. An 
appropriate method is provided by the eikonal expansions, recently introduced in refs. . 
When this procedure is used, a new type of Feynman integrals appears. These integrals and 
methods which were used for their evaluation were described here in some detail. 

The whole construction works well if the mass difference between initial and hnal state 
quarks is not too large. This is the case for the semileptonic 6 —>■ c transitions, where the 
expansion parameter (mb — me)/mb ~ 0.7. In this case we calculated the expansion up to 
the eleventh power of 5; the estimated accuracy of the hnal result is better than 1 %. 

There is, however, a number of other applications where the initial quark is signihcantly 
heavier than the hnal one. It is not obvious to what extend the present method can be useful 
in such situation. There is an indication, however, that our procedures give meaningful 
results even in that limit. As an example, we analyzed the light quark corrections to the 
width of the top quark decay into massless W boson and a 6 -quark. We have shown that 
the hrst several terms of the expansion in St = {mt — mb)/mt approximate the known exact 
result with a 10 % accuracy. 

Finally, we have argued that the same techniques can be applied to the corrections to 
the zero recoil sum rules for the heavy flavor transitions. 
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APPENDIX A: RESULTS FOR THE MAXIMAL RECOIL 


In this appendix we present the hrst nine terms of the expansion of the coefficient func¬ 
tions for the b quark decay rate at maximal recoil, as dehned in (^9]): 


A^ = 
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FIGURES 



FIG. 1. Notation used for the propagators R in diagrams contributing to semileptonic b decays: 
examples with two gluons (a) QED-like and (b) with a nonabelian coupling. 



FIG. 2. The simplest diagram in which the Taylor expansion in {nib — it^c)I'^c leads to spurious 
divergences and must be supplemented by an eikonal expansion. 
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